Lattice effects on the current-voltage characteristics of superconducting arrays 
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The lattice effects on the current-voltage characteristics of two-dimensional arrays of resistively 
shunted Josephson junctions are investigated. The lattice potential energies due to the discrete 
lattice structure are calculated for several geometries and directions of current injection. We com- 
pare the energy barrier for vortex-pair unbinding with the lattice pinning potential, which shows 
that lattice effects are negligible in the low-current limit as well as in the high-current limit. At 
intermediate currents, on the other hand, the lattice potential becomes comparable to the barrier 
height and the lattice effects may be observed in the current-voltage characteristics. 
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Two-dimensional (2D) arrays of weakly coupled 
Josephson junctions in equilibrium are well described 
by the 2D XY model, which exhibits the Berezinskii- 
Kosterlitz-Thouless (BKT) transition driven by the un- 
binding of vortex-antivortex pairs In experiment, on 
the other hand, the systems are usually in nonequilibrium 
and dynamical quantities are measured in the presence of 
external driving [0,1). Although the static nature of the 
BKT transition has been well established, there still re- 
main unsolved questions on the dynamics of the system, 
for example, the value of the exponent of the current- 
voltage {IV) characteristics, the dynamic universality 
class, and the noise spectrum ||^-^. These 2D Josephson 
junction arrays (JJAs) also draw much interest in rela- 
tion with superconducting films and highly anisotropic 
high-Tc superconductors. In JJAs, unlike the latter su- 
perconducting materials, the underlying discrete lattice 
structure causes nonzero pinning potentials for vortices; 
such lattice pinning potentials have not been taken into 
account properly in the theoretical studies based on the 
continuum limit. 

This work investigates the lattice effects in the 2D 
NxN resistively-shunted junction (RSJ) model. We first 
calculate the lattice potential in several geometries in- 
cluding different lattice structures and directions of cur- 
rent injection. The lattice potential barrier on a square 
array with diagonal current injection is found to be much 
larger than that on the same array with horizontal cur- 
rent injection. Comparing the obtained lattice pinning 
potential with the energy barrier which a bound vortex 
pair should overcome to be free vortices, we find that the 
effects of lattice pinning are negligibly small in the low- 
current regime as well as in the high-current regime; this 
is confirmed by the IV characteristics computed numer- 
ically for the square arrays with different directions of 
current injection. It is also found that there exists an in- 
termediate regime of current, where the lattice potential 
effects are observable. 

We begin with the calculation of the lattice potential 
energy, which a vortex should overcome to move to the 



next face of the lattice. In the presence of external cur- 
rents, a vortex is expected to be exerted by the Magnus 
force 

F = -$oJ X z, (1) 

c 

where J is the external current density and $o = hc/2e is 
the flux quantum. Here the direction of the Magnus force 
is perpendicular to that of the current density, enforcing 
the vortex to move along the perpendicular direction to 
the current injection. As an example. Fig. |l| displays 
square arrays in the presence of the external currents 
(a) in the horizontal and (b) in the diagonal direction. 
Under the (sufficiently strong) Magnus force, the vortex 
at position A, which corresponds to the phase configu- 
ration with local minimum energy -Emin7 passes position 
B, with maximum energy i^max, and then moves to C, 
which yields the same configuration as A. Accordingly, 
the potential barrier Eh that the vortex should overcome 
to move is determined by 

I^b -^max -^min (2) 

with the energy 

-E = -Sj^cos(0,-0,), (3) 

where Ej is the Josephson coupling strength, — <j)j is 
the phase difference between sites i and j , and the sum- 
mation is taken over all nearest neighboring pairs. We 
use the method in Ref. Q to obtain the phase configu- 
ration: The minimization of the energy in Eq. (H) leads 
to the condition that the net current fiowing into site i 
should vanish: 

^sin(0,-</.,) = O, (4) 
which can be rewritten as 
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with the summations performed over the four nearest 
neighbors of i. The phase configuration can then be 
found iteratively from 
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where (p^ is the value of at the nth iteration. Equa- 
tion (^), together with appropriate boundary conditions, 
gives the phase configurations for the minimum energy 
£"111111 and for the maximum energy -Emax- In this man- 
ner we compute the lattice potential barrier ioi N x N 
square, triangular, and honeycomb arrays with horizon- 
tal current injection, as well as for a square array with 
diagonal current injection. Figure |^ shows the obtained 
lattice potential barrier Eb versus the inverse size 1/iV, 
which manifests that the barrier heights saturate to con- 
stant values for sufficiently large N. In the case of square 
and triangular arrays with horizontal current injection, 
the barrier heights saturate to the values Eb/ Ej k, 0.199 
and 0.043, reproducing the results obtained in Ref. [||. 
In addition, Fig. ^ also gives the values Eb/Ej « 0.575 
and 0.822 for a honeycomb array with horizontal current 
injection and for a square array with diagonal injection, 
respectively. 

We first consider the energy for unbinding of a vortex 
pair without the lattice pinning potential and then ex- 
amine how the result changes as the pinning effects are 
included. The interaction energy of a vortex-antivortex 
pair separated by distance r is given by EiG'{r), where 
El « 2ttEj at low temperatures and G'(r) is the lattice 
Coulomb Green function (with the diagonal part sub- 
tracted) 0. To an excellent approximation, G"(r) takes 
the form 



G"(r) « \n{r/a)+C 



(7) 



for all r > a, where a is the lattice spacing and C is 
a constant j|] . Accordingly, in the presence of external 
current I the energy of a vortex-antivortex pair reads 0] 



E{r) =Eo + Ei In(^) 
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where Eq is a constant and the last term arises from the 
Magnus force in Eq. (|l[) . We measure the energy in units 
of Ej and write Eq. ^) in the dimensionless form: 



E{r) ^ El lnr-27r/r. 



(9) 



where Eq in Eq. (|^) has been dropped for convenience, 
and r and / are in units of a and the single junction 
critical current Ic = 2eEj/h, respectively. The condi- 
tion for the maximum pair energy dE/dr = gives the 
estimation of the maximum pair size: 



27r/ 



(10) 



for / < £'i/27r w 1 (in the dimensionless form). (For 
/ 1, we have rmax ~ 1-) The energy barrier AE for 
the pair unbinding is thus given by 



AE = E{r 



-E{l)^Ei 
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in the absence of the lattice pinning potential. We display 
AE as a function of / in Fig. ^, where it is observed that 
AE has very large values in the small-current regime. 
Note also that AE = for / > 1, which implies that 
the vortex-antivortex pair can unbind even at zero tem- 
perature if the external current is larger than the critical 
current. 

We now consider the lattice pinning effects on the 
vortex-antivortex pair energy. Taking the position of the 
vortex as the origin, we write the lattice pinning potential 
in the simple form 
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cos 27rr, 
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which is to be included in Eq. (^, with the lattice po- 
tential barrier Eb in units of Ej. In the presence of such 
a lattice potential, the vortex may feel some roughness 
when it moves around. Under small or large external cur- 
rents, however, the lattice effects on the vortex motion 
are expected not to be appreciable: In the latter case of 
large currents, the lattice pinning potential Ep is so small 
compared with the driving potential —2i:Ir in Eq. (|^), 
thus resulting in negligible effects. At small external 
currents, the energy barrier AE for the pair unbinding, 
shown in Fig. || is much larger than the lattice poten- 
tial barrier Eb, and dominates the transport properties 
of vortices since the vortex-antivortex pair should over- 
come the largest energy barrier to be free vortices. On 
the other hand, it is of interest to note that there exists 
the intermediate-current regime, where Eb is comparable 
to AE. In that regime the lattice effects on transport 
properties such as the IV characteristics can presumably 
be measured. 

To investigate the above possibility, we focus on square 
arrays in the presence of external currents in horizon- 
tal and diagonal directions and compute the IV char- 
acteristics at finite temperatures. Since these two cases 
(see Fig. |l|) are believed to differ only in the lattice po- 
tential barrier, we expect that any difference in the IV 
characteristics should be attributed to the lattice pin- 
ning effects. The net current through a Josephson junc- 
tion of shunt resistance R is given by the sum of the 
supercurrent, normal current, and thermal noise cur- 
rent; the resulting current conservation condition at each 
grain yields the equations of motion for a square N x N 
JJA 01 : 
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where If^^ is the external current fed into grain i, the 
primed summation runs over the nearest neighbors of 
grain i, and rjij is the thermal noise current. We here 
employ the fluctuating boundary conditions, and intro- 
duce the twist variables A = (A^,, Aj,) to write the phase 
difference between the nearest-neighboring grains in the 
form 



where is the displacement between i and j and the pe- 
riodic boundary conditions on {(/>i} are imposed in both 
directions. The equations of motion in Eq. ([l^) then take 
the form 01 
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j) + sin(4 - - Tij ■ A) + 



where time has been rescaled in units of h/2eRIc, and the 
thermal noise current in units of satisfies {rjij (t)) = 
and {r]ij{t)VkiiO)) = '2T{S,kSji - 6ii6jk)S{t) with temper- 
ature T in units of Ej/ks- The dynamics of the twist 
variables is governed by two additional equations: 
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k - 4>j - Ax) + riA^ - Ix 
k - ipj ~ Ay) + r/A„ - ly, (15) 
is over all links in the x di- 



where the summation y^/, ;\ 
rection and thermal noise terms satisfy (t?a(0) = 
{VA{t)r]A{0)) = {2T/N^)S{t). In the case of horizontal 
current injection, we have = I and ly ^ 0, while for 
diagonal injection = ly = 1/^/2. By means of the 
Euler algorithm, we integrate Eq. ( |l5| ) and compute the 

voltage V = -y/v^n^, where 14(y) = -N{K^[y))t with 
{■ ■ ■)t denoting the time average. The use of the fluctuat- 
ing twist boundary conditions in the presence of external 
currents has the advantage that the direction of current 
injection can be controlled easily. 

Figure ^ presents the resulting IV curves of the arrays 
of sizes = 4, 8, 16, and 32 under horizontal and diago- 
nal current injections at (a) T = 0.84 and (b) T — 1.30. 
As expected from the existence of the resistive BKT tran- 
sition at T = Tbkt « 0.9, the voltage V atT = 0.84 (be- 
low Tbkt) keeps decreasing as the system size is increased 
in the low-current regime, while V appears to saturate 
to a nonzero value at T = 1.30. It is, however, rather 
difficult to discern the data for horizontal and diagonal 
injection in the logarithmic scale. To manifest the dif- 
ference and to reveal the lattice eflFects in detail, we thus 



plot in the linear scale the difference between the voltage 
under horizontal injection (Vh) and that under diagonal 
injection (V^), which is displayed in Fig. ^ as a function 

(13) of the current / in the system of size = 16 at T = 0.84 
and 1.30. It is obvious that the difference Vh—Vd indeed 
approaches zero in the limit of small and large currents, 
confirming the previous prediction based on the compari- 
son of the energy scales: the lattice pinning potential and 
the energy barrier for pair unbinding. In particular, inde- 
pendence of the voltage in the low-current regime upon 
the direction of current injection suggests that the IV 
exponent for a square array with horizontal current in- 
jection found in Ref. [Q is universal in the sense that it is 
independent of the underlying lattice structure. Further- 
more, since the lattice pinning potential is much smaller 
for horizontal injection {Ei, w 0.199i?,/) than for diagonal 
injection {Ei, w 0.822Ej), vortices with horizontal injec- 
tion can move around more freely, which implies that 
in the intermediate-current regime we have Vh > Vd, as 

= 0, again confirmed in Fig. ^ 

In conclusion, we have studied the lattice effects in 

(14) two-dimensional arrays of resistively shunted Josephson 
junctions for several geometries and directions of current 
injection. The lattice potential energy due to the under- 
lying discrete array structure has been calculated and 
then compared with the energy scale which a vortex- 
antivortex pair should overcome to be unbound. From 
this comparison, we have found that lattice pinning ef- 
fects can be observed in the intermediate-current regime, 
which has been confirmed by the direct numerical inte- 
gration of the current- voltage characteristics. 
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FIG. 1. Illustration of the positions of a vortex under cur- 
rents injected (a) in the x direction and (b) in both x and 
y directions. The positions A and C correspond to the local 
minimum energy -Emin, while B corresponds to the maximum 
-E-max- The dashed linos with arrows describe paths of the 
vortex under the currents. Squares represent superconduct- 
ing islands and arrows denote current injection. 
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FIG. 2. The lattice potential barrier Eb versus the inverse 
system size 1/N . Filled squares, open squares, and filled cir- 
cles correspond to the triangular, square, and honeycomb ar- 
rays with horizontal current injection, respectively; open cir- 
cles to the square array with diagonal current injection. 



FIG. 3. The energy barrier AS for pair unbinding versus 
external current / in the absence of lattice pinning effects. It 
is shown that /S.E = for / > 7c- 
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FIG. 4. Current-voltage characteristics of square arrays of 
sizes AT = 4, 8, 16 and 32 at (a) T = 0.84 and (b) at T = 1.30. 
For comparison, the data for horizontal current injection (la- 
beled by h) and those for diagonal current injection (labeled 
by d) are plotted together. 
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FIG. 5. Difference between the voltage under horizontal in- 
jection (Vh) and that under diagonal injection (Vd) versus cur- 
rent in the system of size Af = 16atr = 0.84 and 1.30. The 
lattice potential effects on the current-voltage characteristics 
are shown to be appreciable only in the intermediate-current 
regime. 
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